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CN ■ Abstract 

3. 

We recast into histories-based form a quantum field theory defined earher 
'nI" I in operator language for a free scalar field on a background causal set. The 

resulting decoherence functional resembles that of the continuum theory 
but the counterpart of the d'Alembertian operator □ is nonlocal and a 
generalized inverse of the discrete retarded Green function. We comment 
bi)! on the significance of this and we also suggest how to include interactions. 
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O ■ 1. Introduction 

O. 

At this stage in the development of fundamental physics one does not know whether 

anything like a quantum field will turn out to exist at the Planck scale, or whether on the 

/S ' other hand quantum field theory will come to be seen as a relatively low energy description 

H : 

of completely different, more fundamental "degrees of freedom" . Once one abandons the 
continuum, a spinor or a vector is no longer a particularly natural object, and in the context 
of causal set theory further difficulties spring from the radical nonlocality inherent in a 
discrete Lorentzian causal structure. Such difficulties suggest that the type of quantum field 
theory embodied in the so called standard model will not survive at a more fundamental 
level. 
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Nevertheless, any attempt to push quantum field theory down to a more fundamental 
level seems bound to be illuminating, and as it happens, we now possess a theory that 
describes a free scalar field propagating in a fixed, background causal set. * On the basis 
of this theory, one can hope to understand better the infiuence of discreteness on such 
phenomena as Hawking radiation or the propagation of wave packets over great distances. 
As formulated so far, however, the theory in question speaks the language of operators 
and state-vectors. The field is described by attaching to each element x of the causet C 
an operator (p{x) acting in a Hilbert space built up from a "vacuum state" by the action 
of creation operators a*. Although all this is very familiar, one might wish also for a 
formulation based on "histories" (a history in this case being a function (/> : C->M). 

Indeed, some of us regard histories-based formulations of quantum theories as more 
basic and more satisfactory than operator formulations, both for the purposes of quantum 
gravity and for the sake of philosophical understanding. How, for example, might one 
expect a dynamical theory of causal sets to look? The only way open to a "law of motion" 
for causets has so far seemed to pass through something like a path integral (or more 
precisely a decoherence functional), and one might say the same about quantum gravity in 
general. Even with the causet limited to the role of "arena" , the knowledge of how to write 
down a path integral for a scalar field on a background causet, would allow the concept 
of "anhomomorphic coevent" to brought to bear on the question of what sort of reality 
quantum field theory is in fact describing [7] . In a more practical vein, one encounters the 
key question of how to generalize from a free field to an interacting one, something which 
so far has not been possible with the operator formulation. As we will see, a path integral 
formulation can suggest a very definite answer to this question. 



For background on causal sets see [1] [2] [3] [4] [5] [6] 

' Sometimes people say "classical history" to emphasize that (pix) is a c-number and not 
an operator. But I will avoid that usage here, in order to avoid any implication that (p 
must obey the classical equations of motion. 



Of equal practical importance is the question of the classical limit. Is it possible in this 
limit to express the equations of motion of a scalar field in such a manner that solutions 
could be built up iteratively given the values of the field in a sufficiently great initial portion 
of the causet? If it were, then computer simulations of wave packet propagation through 
the causet would become possible. On the other hand, if the relevant field equation were to 
involve reference to the distant future, such simulations would become much more difficult, 
or even impossible. In fact an equation of motion of exactly the required, retarded sort 
has been proposed [8] , but in what sense can we claim that it represents the classical limit 
of our discrete quantum field theory? 

In this connection, the following difficulty appears to arise [9]. Suppose that the 
dynamics of the scalar field is to be defined through an action-integral (or sum in the case 
of the causet), and that the integrand at point x depends on the field in the entire past 
of X. When one varies (p{x) in such an expression, the effect will be felt everywhere in 
the future of x. The resulting Euler-Lagrange equations will thus relate (pix) not only to 
the past but to the future. How then could retarded equations of motion result? We will 
return to this issue in sections 4 and 5. 

2. Review of the continuum theory in operator form 

Nothing could be more familiar than a free scalar field in Minkowski spacetime together 
with its vacuum-state and the concomitant Fock representation of the field operators. 
But the way in which these things are ordinarily introduced, starting from the Klein- 
Gordon field equations and the equal-time commutation relations, is not quite suitable for 
transposition to the context of a discontinuous structure like the causal set which not only 
does away with continuous time, but also resists any notion of "Cauchy data at a moment of 
time" . Nor does the normally crucial distinction between positive and negative frequency 



A causal set does admit the notion of slice as maximal antichain E, but such a subset 
unfortunately is not equipped to fill the role of an equal-time hypersurface, or more gener- 
ally of a Cauchy surface in a globally hyperbolic spacetime, the most important difference 
being that the great majority of causal links will "pass thru" E without meeting it at all; 



make much sense in a causal set, which insofar as it resembles a spacetime at all, resembles 
much more a spacetime with curvature than one without it. For such reasons, it has proven 
more fruitful to rest the derivation of the scalar quantum field theory as much as possible 
on the commutation relations alone, or more precisely on the retarded Green function, in 
a way that I will review here. [10] 

By way of preparation, let us first review the theory of a "gaussian" noninteracting 
(massive or massless) real scalar field (f) in fiat spacetime. * 

the usual story 

As this subject is commonly introduced, the input to the theory comprises first of all the 
operator equations of motion in the form of the Klein-Gordon equation, (□ —m )4>{x) = 0. 
Supplementing these equations (plural since there is one for each x) are the canonical com- 
mutation relations, expressed at a fixed moment of time as [(^(x), d(f){y)/dt] = iS^ '{x — y). 
One then expands (p in terms of (four-dimensional) plane waves and identifies the (suit- 
ably normalized) coefficients of the positive-frequency waves with annihilation operators 
a{k). Introducing, finally, the vacuum state |0) via the conditions a{k)\0) = 0, one obtains 
a concrete family of field operators (p{x) acting in a Hilbert space spanned by states of 
the form a* (ki)^^ a* {k2)^^ a* (k^)"'^ ■ ■ ■ |0). From an algebraic standpoint, the equations of 
motion together with the commutation relations define an abstract *-algebra generated 
by symbols (j){x) , while the choice of vacuum induces an irreducible representation of this 
algebra in a particular sort of Hilbert space. The full input to the theory is thus the field 
equations, the commutation relations, and the choice of vacuum (although it is partly 



that is, there will exist pairs of causet elements x, y such that x is to the past of E, y to 
its future, and yet neither E nor any other part of the causet contains an element causally 
intermediate between x and y. In this sense "information can pass through E without reg- 
istering on it" , the impossibility of which is precisely what characterizes a Cauchy surface 
in the continuum. 

with a metric ry^j, of signature ( — h ++) and with \Z\ — v"^ '^a'^b 



a matter of opinion whether or not a change of representation or of vacuum should be 
regarded as changing the theory per se.) 

For our purposes, the most important thing to notice is the need for consistency 
between the field equations and the commutation relations. In the development we have 
just called to mind, this is ensured by the fact that the latter are invariant under the 
time-evolution generated by the former, but we can appreciate better what is going on 
if we express the commutation relations in a manifestly covariant form, the so called 
Peierls form. (The Lorentz invariance of the theory then follows trivially.) To that end, 
let G{x, y) = G^^ {x, y) be the retarded Green function belonging to the wave-operator 
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{n-m')G''\x,y)=S^^\x-y) 



where G is required to vanish unless x y y, meaning that y is within or on the past 
lightcone of x. Further define 

A{x,y) = G{x,y)-G{y,x) , (1) 

which is just the difference between the retarded and advanced Green functions, since (as 
in fact holds in any globally hyperbolic spacetime) the advanced Green function coincides 
with the "transpose" of the retarded one. The commutation relations then assume the 
remarkably simple form, 

[<i>{x),(i){y)] =iA{x,y) . 

If to this equation one applies the wave operator (□ — m ) (acting on x), one trivially 
obtains zero on the left hand side. The required consistency then follows from the fact 
that the right hand side also vanishes because A is the difference of two Green functions 
for (□ —m ). Notice here that the last three equations all make perfect sense, and are 
valid, in any globally hyperbolic spacetime. The resulting definition of a quantum field 
theory accordingly goes through as well, except that something is needed to replace the 
notion of positive frequency if one wishes to define a vacuum state. 



Now that we have a field theory in operator form, the next question is how to render it 
into histories form. For this there is a standard procedure that I won't review here, because 
we wiU meet it again in detail when we come to the case of the causal set. Instead, let 
me just quote the resulting decoherence functional, which by definition attaches a complex 
amplitude to every pair of spacetime histories of the scalar field. (A pair of histories can 
be called a "Schwinger history" in honor of the so called Schwinger-Kel'dysh version of the 
path integral, which operates precisely with such pairs.) To avert confusion between (p as 
a ("classical") history and (p as an operator- valued field, let me for now use the symbols ^ 
and ^ to represent histories. The decoherence functional D{^,^) is then given formally by 
the expression 

exp{i5(0} exp{z5(0}* 5(e|Sr,?|Sr) (2) 

where S{^) = \ /( — (V^) —m ^ ) plus an implicit further term that incorporates the effect 
of the "initial state" , in this case the vacuum. Also E j- is a spacelike future boundary at 
which the spacetime has been truncated and the 5-function involving it has the effect of 
forcing the otherwise independent histories ^ and ^ to share the same restriction to this 
final boundary. From this, the decoherence functional extended to sets of histories and the 
corresponding quantum measure can be obtained in the usual way [11]. 

For future reference let us also define here the so called Wightman or "two-point" 
function given by the expectation of (/){x)(/){y) operators in the vacuum state: 

Wix,y) = {<P{x)<Piy)) (3) 

and let us note also that in the vacuum 

W{x) := {(p{x)) = . (4) 

Because the vacuum is a "gaussian state" the entire field theory is effectively encapsulated 
in these two equations, a fact we will use heavily in our construction of the decoherence 
functional in the causet case. 



a different starting point 

The logic of the "usual story" that we have just rehearsed is simplest if we present the 
commutation relations in their Peierls form. We might then summarize in part the steps 
we followed as: wave operator — )■ G^^ — t- A — )• [^, ^] . But consistency between the wave- 
operator, □ —m , and the consequent commutation relations required, as we saw, that 
the transpose of the retarded Green function also be a Green function for the same wave 
operator. Since no such relation is known in the causet case, it turns out that one can 
get farther by changing the starting point slightly from the equations of motion to the 
retarded Green function G^^ . From it we can derive the commutators in just the same 
way. Indeed, one can pass freely from G^^ to A and back by means of equation (1) and 
the inverse relation G^^ (x^y) = A(x,y), which holds when x y y. (In equation form, 
G^^^{x, y) = 9{x, y)A(x, y) where 9 is the "covariant Heaviside function".) 

But how to continue without the aid of any equations of motion? The surprising 
answer discovered implicitly in [10] is that one can pass directly from the commutator 
function A to the two-point function W{x, y) by thinking of A(x, y) as a real skew matrix 
A^^ and forming the positive part of the corresponding Hermitian operator zA. Our new 
logic will therefore go as follows: G^^ — )■ A — )■ W . It would be difficult to imagine anything 
much simpler! 

Notice in this connection that W always is, in fact, a positive (semidefinite) operator, 
as follows immediately from the fact that v*Wv > for any complex vector v. Indeed, 
positivity of the Hilbert space inner product yields for any such v^ , '^{v^)*W^^Vy = 
Y.{v^y{(t)'"(t)^)Vy = {X^X) = \\X\0)\f > , where X = E^x'/'''- It's equally true that W, 
the complex conjugate matrix of W, is positive semidefinite and that iA = W — W, as 



follows from: {W - W)""" = (0"(/)^) - {(^^(py) = ((/)"(/)^) - ((/)^(/)") = ([(/)", (/>^]) = (zA"^) = 
zA^^ Thus, zA is always the difference of two positive matrices W and W. What our 
prescription for W is adding to this fact is that these operators have orthogonal support: 

WW = WW = . (5) 



One might view this as a kind of "ground-state condition" imposed on W beyond what 
foUows automaticaUy from the fact that it is the two-point function of a selfadjoint operator 
(p^ . From (5) follows immediately a simple matrix equation giving W in terms of A (and 
therefore in terms of G) : 



M^ = l/2(zA+ V-A2) . (6) 



In practice one would often compute y/—A'^ by diagonalizing A. Once this was done, 
one would obtain W just by expanding iA in terms of its eigenvectors and retaining the 
terms with positive eigenvalues, and one could regard this as the practical meaning of our 
prescription for obtaining W from A. 

For future reference let me introduce the notation R = ^ V— A^, a manifestly real and 
positive matrix. Comparison with (6) then shows that R is just the real part of W, whose 
decomposition into real and imaginary parts is thus given by 

W ^R + iA/2 . (7) 

It's worth noticing here that the definitions we have given yield a free field theory 
with a distinguished "ground state" for any globally hyperbolic spacetime or region of 
spacetime. (Modulo, however, certain technical questions about convergence that arise 
when the region's volume is infinite.) Notice also that nowhere do these definitions refer to 
any notion of positive frequency. In Minkowski spacetime, the new prescription can easily 
be seen to reproduce the usual vacuum (it more or less had to by Lorentz invariance). In 
this sense it generalizes the Minkowski vacuum to the case of arbitrary curvature. 

Perhaps also it's also worth mentioning a possible further generalization that would 
tie W less closely to G. In order to reproduce correctly the commutation relations, it is 
necessary only that the imaginary part of W coincide with A/2. It is also necessary for 
mathematical consistency that W be positive, and perhaps other conditions are needed as 
well, but it is clearly not necessary that W simply be the positive part of iA. One could 
thus consider relaxing condition (5) as well as the assumption that {(p) = 0. As will be 



seen, neither generalization would require major changes to our discussion below of the 
decoherence functional on a causal set. 

3. Scalar field theory on a causet in operator form 

Let us now leave the continuum behind and turn to the case of a causal set. Henceforth our 
discussion will be in the context of a fixed, finite causal set C, and it will assume further 
that a fixed "retarded Green function" has been adopted, which I'll write as either G{x,y) 
or G , where x and y vary over the elements of C. In some of the most important cases a 
convincing candidate for G is known and has been tested in practice to a greater or lesser 
extent. These cases include zero mass in a causet intended to approximate an arbitrary 2- 
dimensional spacetime, as well as arbitrary mass in a causet intended to approximate a flat 
spacetime of dimension either 2 or 4. (See [12].) In addition, the "retarded d'Alembertians" 
of [8] and [13] could be inverted to produce candidate Green functions in more general 
causets, but this type of prescription has not yet been put to a test. 

Starting from G and following the same steps as delineated above yields then a scalar 
field theory for the given causet, complete with operators (j){x), vacuum-state |0), and 
Wightman function W{Xjy) = {(/){x)(/){y)) , where all expectations are taken with respect 
to the state |0). 

As we have already noted, the fact that this theory-cum-state is in some sense gaussian 
means that all its consequences can be recovered starting from W alone. It might be 
interesting to tease out better in exactly what sense the word "gaussian" can be applied to 
our causet theory. Having no Cauchy surfaces, it also has no straightforward Schrodinger 
representation, and therefore no vacuum wave function whose exponential form could be 
held up as the meaning of the word. On the other hand, the fact that the (p^ are linear 
combinations of raising and lowering operators implies that Wick's theorem will work in 
the expected manner, and this is all that we will need. Indeed, we will need only the 



following consequence of Wick's theorem, where $ stands for any linear combination of 
the field operators (j){x). 

(exp{$})=exp|i^| (8) 

In view of equation (20) of the appendix, the following calculation demonstrates this 
identity. 

n n n ' n 

This seems a good place to mention a way of thinking about W that is more geometri- 
cal in nature than the characterizations we have encountered so far and that involves only 
real, as opposed to complex, vector spaces. Namely W can be represented geometrically as 
a family of mutually orthogonal, "axially weighted 2-planes" in the real vector space M , 
where A^ is the cardinality of the underlying causet, these being a family of pairs {w, A) 
where w is an oriented 2-plane in R and A is a strictly positive real number corresponding 
to one of the eigenvalues of iA. This one can see by splitting the eigenvectors of iA into 
their real and imaginary parts, or equivalently by decomposing A as a sum of terms a /\h 
with a^h E ^ being orthogonal vectors. The latter vectors then provide the basis of a 
singular-value decomposition of the real matrix A. Thanks to these relationships, one can 
construct W directly from A by forming the singular value decomposition of the latter, 
thereby economizing — in computer simulations — on memory and CPU-time, since only 
real numbers need be involved. 

We are almost ready to turn our efforts to producing a histories version of our theory, 
but first, perhaps a simple example of the scheme G — )■ A — )■ VF would be in order. Consider 
then a causet of only two elements, e^ -< e-^, making up a 2-chain. (Given the interpretation 
of VF as a collection of weighted 2-planes, we could in this very simple case pass immediately 
to a unique W without bothering with G and A, but that would not illustrate the general 
situation.) That G is retarded means precisely that G =0, but since only off-diagonal 
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terms survive the antisymmetrization in (1), we might as well suppose that the diagonal 
of G vanishes as well. Up to sign and normalization, G is then unique: 

fO 
^-^1 

from which there result immediately A = j J and 

(the Pauli matrix). Substituting this into (6) then furnishes W as 



l^^-^^-Ul T 



W = -(1 + (T2 



i.e. 2R in (7) is simply the identity matrix. We could also have arrived at the same 
result by diagonalizing iA explicitly as cr2 = ^^ ~ uu\ where u = l/v/2 I • I , and then 
discarding the negative term. Notice finally that W correctly satisfies the "ground-state 
condition" (5). Given that W is hermitian, (5) is just the assertion that WW = 0, tilde 
denoting the matrix transpose. Alternatively expressed, this says that the dot-product 
(the bilinear one, not the hermitian one) of any row (resp. column) with itself or with any 
other row must vanish, and this is true by inspection. 

4. Scalar field theory on a causet in histories form 

Henceforth it will be convenient to redefine G = G-' to be the retarded Green function 
with its diagonal set to zero. Since A is G minus its transpose, it is unaffected by this 
modification, and since W was defined through A, it is not affected either. Equivalently, 
we could take W as our starting point and then simply define G to be twice the retarded 
half of the imaginary part of W. Either way its clear that for purposes of this paper, there's 
no harm in dropping diag(G). But in a larger context this change would need to be borne 
in mind. For example, if we had obtained G by inverting some retarded D'alembertian 
operator for the causet, then we could not recover the latter from G without restoring its 
diagonal. 

11 



By a histories-based formulation of a field theory, one means (for a real scalar field) a 
formulation that works directly with field configurations (p : C ^W { "histories" ) , and avoids 
reference to field operators or state-vectors, except possibly as auxiliary technical devices 
[14]. In practical terms, this means a path-integral formulation, although one might hope 
that more general frameworks will one day be available. In order that a framework based 
on path-integrals be self-contained, moreover, it seems necessary to assign amplitudes to 
Schwinger histories rather than to individual histories, as the more familiar path-integral 
does. That is to say that one needs to express the dynamics in terms of a decoherence 
functional D, from which the corresponding quantal measure [15] can then be computed. 
Since the definitions of these things have been presented many times in the literature, I'll 
not repeat them here. Rather, I'll begin by writing down the expression (an expectation 
of a product of projectors) that one needs to evaluate in order to recover the decoherence 
functional from a theory expressed in the operator language. 

In line with the definitions laid down above, let ^ and ^ be two (completely indepen- 
dent) histories, each being specified by a list of real numbers ^^, one for each element x of 
the causet. The complex number D{^,^) is then given by the equation 

Di^. = {^i^' - tn<p' - f ) ■ ■ ■ sicp"" - fncp"" - e^) ■ ■ ■ 5{<p' - en<p' - e)) m 

Here (/)■' = (j){xj), ^■' = Ci^j)^ etc, where j = 1---A^ is any natural labeling of C. In 
other words the elements x G C of the causet must be labeled so that no element with a 
smaller label ever temporally follows an element with a higher label: x -< Xf. ^ j < k. 
Such a labeling always exists, and any two natural labelings are guaranteed to produce 
the same result for -D(^, ^) because A(x, y), and hence [4>{x), </>(?/)], automatically vanishes 
when X and y are causally unrelated ("spacelike"). Notice that D is normalized such that 



perhaps obtained by freeing the concept of preclusion from its dependence on the no- 
tion of quantum measure and measure-zero, or perhaps even by freeing anhomomorphic 
inference from its dependence on preclusion, cf. [7] 
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Re-expressing the ^-functions as integrals puts (9) into the form 



which represents -D(^, ^) as the Fourier transform of what one might caU its "non-commuta- 
tive characteristic function" in analogy with the concept of characteristic function that 
figures in ordinary probabihty theory. (Here again the real numbers A and A • are inde- 
pendent parameters of the Fourier transform. The bar is not being used to denote complex 
conjugation.) 

Before plunging into the evaluation of this expression for the decoherence functional, 
let us evaluate the simpler expression which results from integrating out the ^ and all but 
one of the variables ^ = ^(x) in -D(^, ^), which mathematically is analogous to a marginal 
probability density for the remaining variable (f){x). After doing the integrals, we are left 
simply with 

— exp{-zAO(exp{zA(/)(x)}) , 

which in light of the identity (8) turns into 



|^exp{-«;V/2-zAe} 



y2Tiw 

where I have written ^ for ^{x) and w for {(p(x)(p{x)). Unsurprisingly for a free field, we 
just obtain a gaussian. 

The next simplest case, which already illustrates most of the complications, results 
from integrating out all but two of the field- values, and yields an integral of the form 

/ 7^7^exp{-za^-i/377}(exp{za(/)(x)}exp{z/3(/)(|/)}) . 
J Ztv Ztv 

We can again avail ourselves of (8), but only after bringing to bear the well known identity 

^A^B^^A+B^l[A,B] ^ (11) 
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which holds whenever the commutator [A, B] is a c-number. In virtue of (11) we find 

exp {ia(l){x)} exp {i(3(f){y)} = exp < ia(l){x) + ij3(l){y) — -a(3[(l){x), (l){ii)\ 

which can also be written as exp \ia(l){x) + i/3(/){y) — ^al3{W^^ — VF^^)} because 

mx),<Piy)] = <Pix)<p{y) - (PiyMx) = {<P{x)<p{y) - <Piy)<P{x)) = W^^ - W^^ 

the original commutator being a c-number. Hence, 
(exp {ia(/){x) } exp {z/30(y) }) 
- exp |-1 {{{ac^ix) + my)f) + «/5(W^"" - W^"")) 

= exp |-^ (aV"" + (3^Wyy + afSW^ + a(5W'''" + ^(5^^ - a(5Wy'') \ 

= exp |-- (a'W^"" + p'^Wyy + 2a/3Ti^"^) | 

= exp|-i((a p){^^y wyy){li))] 

The important thing to notice here is that the order of the indices in W^^ copies the order 
of X and y in the original factors. 

Following this pattern, it is straightforward to evaluate the quantity within angle 
brackets in (10). The result is 

exp |-1 (#^%A, + W^'^W + 2W^%\) I (12) 

where "^ ^^ denotes W^^ li j < k and W^^ li k < j, for example #^3 ^ ^32 ^ ^23. 
and W follows the reversed convention. (That the arrows in the first and second terms 
point in opposite directions merely refiects the opposite ordering of the causet elements in 
the two halves of the Schwinger history.) Notice in connection with (12) that the matrices 
W and W are complex conjugates of each other, and recall also that W itself is positive 
semidefinite. 
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Let us observe in passing that integrating out the ^ variables from D{^,^) would have 
the effect of setting A = in (12), leaving in the exponent only —^W-' XjXj., which involves 
only the time-ordered Wightman function, or equivalently the "Feynman Green function" . 

Now since D{^,^) is nothing more than the Fourier transform of (12), its evaluation 
would be relatively carefree, were it not for the fact that the quadratic form in (12) is not 
invertible. This will prevent D{^,^) from assuming the pure gaussian form that it would 
otherwise have had, and that one might naively have expected from a free field theory. 
Instead, it will be the product of a gaussian with a 5-function that will enforce certain 
"constraints" on its arguments. 

The further manipulations to be done will become more transparent if at this point 
we switch to a matrix notation, and I'll continue to use a tilde to denote transpose. Then 
for example, we can re-express (12) as exp { — 1<5} where 

Q = XW X + XV^X + 2XWX . 

Let us now express Q directly in terms of G and R using that, according to its definition, 
W = R + zA/2 = R+ {i/2){G - G), and noting also that ^ = ^ =- i? and ^ = -^ = 
G + G, whence 

W = k + iA/2 = R + iA/2 = R+ {i/2){G + G) 

whence similarly W = R— {i/2)(G + G). Expanding W this way yields 

XW X = X{R+ {i/2){G + G))X = XRX + iXGX 

with similar expansions for XW X and 2XWX. Adding these together then reveals that Q 
is expressed most simply in terms of sum and difference variables for A and A. 

g = (A + X)R{X + A) + z(A + A)G'(A - A) . 

Its fourier transform can then be expressed in terms of the corresponding sum and difference 
variables for ^ and ^. To that end, let us define 

K = X + X,k = ^{X-X), </,= l(^ + ^), y, = ^-^ 

15 



Here (/>, the mean of the two halves of the Schwinger history, is what is sometimes caUed 
"the classical field" because ^ = ^ in the classical limit, while v' represents in some sense 
the deviation from classicality. 

With these definitions, A^ + A^ = Kcj) + A;(/7, and the fourier transform yielding -D(^, ^) 
turns into 

I d^K d^k exp I -iKcP - ikip - -KRK - iKGk \ (13) 

where I've left out an overall normalization which needn't be carried along because it can 
be recovered at the end from the condition f d ^d ^D{^,^) = D{Q, O) — 1. 

Evidently the conjugate variable k occurs only linearly in the exponent of (13), in the 
combination k(f + KGk = k{GK + </?); the dk integral therefore produces the 5-function 
I mentioned earlier, which constrains </? to be in the image of G, the "advanced Green 
function". Equivalently, it constrains </? to be orthogonal to ker(G) = {K \ GK = 0} . We 
now have 

D{^,0 = {const) f d^K S{GK + ip)expl-^KRK-iK(j)\ . (14) 

When (f is in the image of G, as required for D to be nonzero, (14) has the form 
of a gaussian integral over an affine subspace of M , namely the space of solutions to 
the equation GK + ip = 0. To convert (14) into an integral over a full vector space, let 
/ : y — 7- R be the inclusion map of V^ = ker G into M and let Kq be any solution of 
the equation GKq = — </? , so that the space over which we have to integrate is Kq + V. 
(Concretely we can represent / by a matrix / ", where the vectors /•" for fixed a furnish 



Please don't confuse (p as defined here with the same symbol used earlier for the field 
operator. 
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a basis for V.) The substitution of Kq + K for K then yields our integral in the form 

/ d^' Ke^Y>\--KRK-iK(j)\ 

JkeKo+v 12 J 

= / d^'Ke^p\-hK + Ko)R{K + Ko)-i{K + Ko)(j)\ 
JKev I ^ J 

= exp I -]-KqRKq - iK^A f d^' K exp | -]-KRK - K{RKq + i<P) \ 
where A^ = dimkerG is the dimensionality of V. Then with the aid of the matrix / we 

N' 

obtain an integral over M by putting Iv in place of K : 

exp I --KqRKq - iKqCJ) I / d^'v exp | --vIRIv - {K^R + i(l))Iv \ . (15) 

The result is now an ordinary gaussian integral over M of the form J exp \^ — ^vAv + Bv^ 
with A = IRI and B = {KqR + i(p)I. Under the assumption that A is invertible, this 
integral evaluates to (a constant multiple of) exp | — |i?(l/A)i?} . When expanded out, 
the exponent here consists of three terms: 

(j)I{IRI)~^I(l), (16-1) 

-iKo{l-RI{IRI)~^I)(f) , (16-2) 

--Kq {R - RI{IRI)-^IR) Kq . (16 - 3) 

We now have the decoherence functional in hand, except that (16) is expressed in terms 
of the ambiguous constant Kq rather than </?. For consistency, this apparent dependence 
on Kq must be illusory, i.e. (16) must be independent of Kq for fixed </?. 

In considering (16), let us observe first of all that one would expect / to be "small" 
since one would expect dim ker G to be small compared to A^, the number of elements in our 
causet. It can't vanish entirely, given that G is strictly retarded — now that we've removed 
its diagonal — but an analogy with the continuum might suggest that it would receive con- 
tributions only from the "boundary region" of the causet. To the extent that / actually can 
be neglected, (14) could have been evaluated without doing any integrations at all, as the 
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5-function therein sets K to —G~ </?, yielding D = exp < —^ip{G~ RG~ )(p + i(pG~ (f)>. 
Now a natural symbol to represent G~ would be D 5 since in the continuum G is in some 
sense the inverse of □ . In the present setting we can also appreciate that G~ (if it actually 
existed) would be retarded too because G itself is retarded. One might then refer to □ = 
G~ as a "retarded d'Alembertian" and its transpose □ as an "advanced d'Alembertian" . 
With this nomenclature we would have D = exp < —^^{\Z1 Rlll)^ + '>'^ \Z\ (p} ■ Moreover, 
in the continuum, \Z\ R = Re(n W) = 0, therefore one might also expect —^(p{\Z\ Rlll)^ 
to be small in the causal set. In that case the whole of the decoherence functional would 
simplify to exp {i(f □ (p} . Remarkably, this is exactly the decoherence functional of the 
continuum theory, as an integration by parts reveals, given that </? = on the future bound- 
ary. More precisely, the continuum decoherence functional is this expression together with 
the 5-function in (2) that mandates the just-mentioned equality between ^ and ^ on the 
future boundary, cf. [11]. Such a 5-function is present implicitly in (14) and explicitly in 
(9), but I don't know how to extract it as such from a well-defined continuum limit of the 
discrete theory. We'll return to the constraints on </? shortly. 

It is useful at this point to distinguish conceptually between the space M of field- 
configurations and the space M of configurations of the dual variables K, A, etc. Since 
we have already employed the symbol O for the former, we might as well denote the latter 
by O* , using a star to denote dual space. I'll also write V for R considered as kerG 
considered as a separate vector space. Then I : V -i-Q* , I : Q^V* , G:0*— t-O, 
5:0*^0, etc. 

Let us return now to the full decoherence functional as given by (16). Its second 
member can be written as —iKQ{l—Il)(j) if we define 11 to be the operator 11 = RI{IRI)~ I . 
Here, of course, we continue to assume that {IRI)~ exists. From its definition it is trivial 
to check that H, though not hermitian, is a projection in the sense that 11 = H; and also 
that 

7n = 7, UI = I. (17) 
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Now the term by which Kq = —G~ cp is ambiguous is by definition Iv for arbitrary 
V E V. Hence the ambiguity in —Kq{1 — lV)(p is precisely {Iv){l — Ii)(p = vl{l — Il)(p = 
v{I — IIl)(p, which vanishes by (17). Consequently, it is justified to write (16-2) simply as 
i{G~ f){l — n)(^ = iipG~ (1 — n)(/) . In view of these considerations, it is natural to define 

n = G-'(i-n), (18) 

meaning that □ solves the equation G \Z\ = 1 — H. We have just seen that □ (/> is 
unambiguous when contracted with any ip G im(G). Thus □ itself is uniquely defined if 
we construe it as a map from O to (iniG)*, the dual vectorspace to (imG) : in symbols, 
□ : O — 7'(imG)* C O* . Regarded however as mapping O to all of O*, □ is ambiguous by 
the addition of any linear operator with image in (im G) = ker G . A nice idea at first sight 
would be to take advantage of this ambiguity to render □ fully retarded, but unfortunately 
that is not possible. Reasoning similar to the above also proves that G \^ G = G, that 
imll = im/ = kerG, that ker 11 = imG, and therefore that im(l — H) = imG. Also that 
UR^RU = RI{IRI)~^IR . 

So far, we have been occupied primarily with the analysis of (16-2), but we also need 
to deal briefiy with (16-1) and (16-3). The former is complete as it stands, but the latter 
requires to be expressed in terms of (/?, just as we did for (16-2). To that end, rewrite (16-3) 
as -^Kq{1 - U)RKf^ = -ii^o(l - ^^RKq = -|^o(l - n)i?(l - U)Kq and observe that 
the ambiguity of adding a vector Iv to Kq drops out just as before because (1 — H)/ = 0. 
Continuing on, and replacing Kq by —G~ </?, we obtain with the aid of equation (18), 

— ^ipG~ (1 — n)i?(l — n)^" if = —\^ n(l ~ n)i?G(/7, which can also be written as 

— \{\I\^)R{\I\^)- Collecting our three terms then yields in total 

D{i,l) = {constant) exp jiv? D </* - ^ (/</>) (7i?/) "'(/(/)) - \{Uv)R{nv)\ (19) 

In view of the rather intricate derivation we've just been through it seems worth 
recording here the much simpler expression that results from integrating out the ^ variables 
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in -D(^, ^) or simply from returning to (9) and deleting the factors involving ^. One sees 
immediately from (12) that — provided that W is invertible — the result will be 

D{C,n) = Jd^'^DiU) = (const) exp l^e D ^1 

where I've written Q^ for iW ~ , a notation that seems natural because —iW corresponds 
to the Feynman Green function of the continuum theory. 



5. Consequences and questions 

Perhaps the most important conclusion from the above work is that the equations of mo- 
tion for the "classical field" (p are essentially a retarded non-local version of the continuum 
equations □ (/^ = 0. Equations of precisely this sort were proposed in [8], but until now 
there has seemed to be no sound reason for preferring retarded equations over, for exam- 
ple, the time-reverse or some combination of the two. Granted that retardation greatly 
facilitates computer solution, but that would be more an opportunistic reason than one 
of principle. Retardation also agrees better with the notion that physics in the classical 
regime, even if it were to be non-local, would still be "causal" , but here we have reached 
a similar conclusion in a far more convincing manner! 

But why should one refer to □ (^ = as the classical equations of motion? Expressed 
through the decoherence functional, classicality amounts to the vanishing of ip, i.e. to the 
equality of ^ with ^ (in which case both coincide with (p). But if we vary ip in (19) and then 
set (/? = 0, the result is precisely □ (/> = 0. Further support for this point of view comes 
from the fact, already alluded to, that the second and third terms in (19) seem likely to be 
relatively negligible, having no further significance than either initial-conditions or small 
corrections. In the case of the second term specifically (the term in (p) , its character as a 
kind of initial condition or "state" can be perceived fairly clearly, because it depends on 
(p only through I(p, which is a combination of contractions of (p with vectors in ker G. But 
since G is "future-looking" it will annihilate any function with support on the minimal 
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level of C. Plausibly then, its entire kernel is comprised of functions supported near this 
bottom level, whence Icj) would depend only on such initial values of (j). 

Of course if we vary all the arguments in (19), we obtain more equations, and taken 
together they tell us more than just □ (/> = 0. In fact they constrain (j) itself to be zero 
(along with </?). This however only reflects the fact that our action (or rather decoherence 
functional) implicitly includes a specification of "initial conditions", something like an 
"initial state wave-function" in the continuum theory. In fact exactly the same thing 
happens in that case: a unique "classical" solution is picked out. In the present context 
it's unremarkable that the corresponding solution is (/> = 0, because our construction of 
the decoherence functional assumed early on that (4) was true. 

We might also ask how an apparently time-symmetric scheme has given rise to a 
very time-asymmetric set of classical equations. To the extent that we can regard □ 
as retarded, the combination «</? □ (^ is entirely retarded in the "classical component" cj) 
and advanced in the "non-classical component" </?. Just how this has come about seems 
mysterious, but logically the asymmetry had to originate in the order chosen for the factors 
in (9). Perhaps a careful tracing through of the steps leading from there to (19) would 
shed more light on the question. 

Of course the operator □ is not fully retarded but only approximately so (if our 
plausibility arguments are valid). Should this worry us? I'm not sure, but one reason for 
equanimity at this point springs from the observation that a dynamics taking place in a 
frozen causal set can only be a partial refiection of full quantum gravity. Really the causet 
itself needs to be included in the decoherence functional. Not having done this we have 
restricted its future development in an unphysical manner, and thereby we have imposed 
a sort of future boundary condition. It would be no surprise if such a condition were to 
be refiected to some extent in the dynamics of the variables (namely the scalar field) that 
we have left free. (Whether we can be equally happy with an anti-causal equation for tp is 
another question however.) 
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A question might be raised here about our assumption that the matrix IRI is invert- 
ible. How well does this hold up in practice? Quite well if the few numerical experiments 
that I've done are indicative, but so far they have been limited in scope. An interesting 
exception, however, springs from the existence of so called non- Hegelian elements or pairs 
within the causal set C, viz. pairs of incomparable elements x and y which share the 
same relation to every third element z E C. Such elements cannot be distinguished from 
each other by the causal structure, and in consequence the "difference function" defined 
by f{x) = +1, f{y) = —1, and f{z) = otherwise, will be in the kernel of G, A, etc. 
In particular it will be in the kernel of IRI, rendering the latter non-invertible. But the 
problem this might seem to raise solves itself. Indeed x and y are in some sense not two 
distinct elements at all, and the theory seems to know this, as one can see by returning to 
equation (13) or (14) and noticing that / is a degenerate direction for the quadratic form 
involved. Consequently the resulting decoherence functional will vanish unless ^{x) = ^{y) 
and ^(x) = C{y)- The theory thus lives in effect on the quotient causet C' obtained by 
identifying x with y, and one might just as well eliminate y (say) at the outset and continue 
in this manner until no further non-Hegelian elements remain. 

Many more questions could be raised at this point, but let me just conclude with 
two of a more general nature, the first being this. All our work herein has presupposed 
a fixed, background causal set C with a finite number of elements. But in reality, one 
expects the causet to be growing and infinite, or at least potentially infinite toward the 
future. We cannot actually handle dynamical growth without setting up a full quantum 
gravity theory, but we can certainly imagine C being infinite toward the future. In that 
case our theory herein would correspond to only some initial portion of C (a so called 
partial stem) and this partial theory should be obtained from the full one by restriction. 
Conversely, if the theory is formulated initially on the stems, then the need arises for 
consistency between the partial theories belonging to the different stems of C. (From a 
histories point of view, almost all quantum theories are built up this way from partial 
theories which provide the measures of so-called cylinder events. See for example [16] or 
[17].) It is therefore an important question whether the theory herein can be, or can be 



22 



adapted to be, a component of a coherent theory defined consistently on an untruncated 
causal set. 

The other general question concerns interactions. The theory we have been working 
with herein is that of the causet counterpart of a free field, as manifested in the quadratic 
nature of (19). In its original operator form, however, the theory does not readily suggest a 
generalization to include interactions. A potential advantage of the histories formulation, 
then, is that it does suggest such a generalization. To incorporate a (p interaction, for 
example, one need only include a multiple of i^ — i^ into the exponent of (19). For 
the sake of consistency, one should check that the decoherence functional, thus modified, 
remains positive semidefinite, but this can be done. It would be interesting to compare 
the consequences of this generalized theory with those of its continuum counterpart. 

Appendix. Proof of an identity used in the text 

In proving equation (8) we used the following identity (20), which is a special case of 
Wick's theorem. Indeed, the coefficient {n — 1)!! = (n — l){n — 3) ■ ■ ■ 5 ■ 3 ■ 1 in (20) is 
easily seen to count the number of ways of pairing up the 2n factors of $ in the expression 
($"^) = ($$$ ■■■$). Here we prove this special case directly by a straightforward method. 

($2-) = (2n-l)!!($$)" (20) 

Suppose first that ^ = a + a* with a normalized as usual so that [a, a*] = 1. Also let |0) 
be the "vacuum" relative to a with (■) = (0| ■ |0). Noticing that a|0) = and (0|a* = 0, 
and that [a, $] = [a,a + a*] = [a, a*] = 1, we have then 



($-) = ((a + a*)$"-^) = (a$"-^) = ([a, $"-^]) = {{n - 1)$"-^) = {n - l)($"-2) 



where we used the fact that [a, •] is a derivation (it obeys the Leibniz rule). Since ($ ) = 1, 
it follows by induction that for even n, ($'^) = (n — 1)!!. When n is odd ($"^) of course 
vanishes, since an odd number of applications of a or a* to the vacuum can never bring 
one back to the vacuum. 
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Finally, consider the more general situation where $ is some linear combination of 
field operators (p{x) . Because $ = $* we can always write it in the form (p = A + A* , 
where A contains only the terms with lowering operators a. The proof then proceeds as 
before except that in place of [a, a*] = 1 we have instead [A, A*] = ($$), as we see from 
($$) = {{A + A*){A + A*)) = (AA*) = {[A, A*]) = [A, A*] . 
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